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Chemical defects in polyethylene (PE) can deleteriously downgrade its electrical properties and performance. 
Although these defects usually leave spectroscopic signatures in terms of characteristic luminescence peaks, 
it is nontrivial to make unambiguous assignments of the peaks to specific defect types. In this work, we go 
beyond traditional density functional theory calculations to determine defect-derived emission and absorption 
energies in PE. In particular, we characterize PE defect levels in terms of thermodynamic and adiabatic 
charge transition levels that involve total energy calculations of neutral and charged defects. Calculations 
are performed at several levels of theory including those involving (semi)local and hybrid electron exchange- 
correlation functionals, and many-body perturbation theory. With these critical elements, the calculated 
defect transition levels are in excellent correspondence to observed luminescence spectra of PE, thus clarifying 
and confirming the origins of the observed peaks. Based on this work, a prescription with a reasonable 
computational expense is proposed to accurately predict and assign spectroscopic signatures of defects in 
other organic polymers as well. 


I. INTRODUCTION 

Polyethylene (PE) is an important insulation material 
that has found widespread use in electrical applications 
including transmission line cables and capacitorsf 1 The 
electrical performance of PE over the long term is af¬ 
fected by impurities and chemical defects that are origi¬ 
nally part of the material, as well as those that are cre¬ 
ated progressively with time. Such defects can introduce 
charge carrier (or defect) states within the band gap of 
PE, can act as traps and sources of charge carriers, cat¬ 
alyze further damage, and can deleteriously affect the 
overall conduction behavior of the insulator ! 2 ^ It is thus 
critical that a firm understanding of the nature of such 
defects be obtained. 

The best evidence for the presence of chemical de¬ 
fects (especially in PE), and a knowledge of their type, 
is provided by a variety of luminescence measurements, 
including electroluminescence, photoluminescence and 
chemiluminescence P ^ 9 Defects, depending on their type 
and the specific details of the energetic placement of their 
defect levels, lead to characteristic emission signatures. 
Such measurements have led to the identification of a 
plethora of defects in PE, the chief among them being 
the carbonyl, dienone, vinyl, hydroxyl, etc., as illustrated 
in Fig. [l] Nevertheless, assignments of the luminescence 
emission bands are never straightforward, and an alter¬ 
nate check of whether a particular defect will lead to a 
particular emission peak or band is highly desired. 

Such an alternate route to identify emission signatures 
can be provided by first principles computations, e.g., 
those based on density functional theory (DFT )® 11 In¬ 
deed, DFT and beyond-DFT computations have played 
a critical role in successfully unraveling the repercus¬ 
sions and signatures of defects in inorganic semiconduc- 
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FIG. 1. Six typical defects that occur in polyethylene are con¬ 
sidered in this work, e.g., carbonyl, dienone, hydroxyl, double 
bond, conjugated double bond and vinyl. The red, white and 
grey spheres are O, H and C atoms, respectively. 

tors and insulators in the past JEHH] i n the case of PE, 
such work is in a state of infancy .® 19 Past DFT work 
on PE is dominated by the use of (semi)local exchange- 
correlation functionals to treat the quantum mechanical 
part of the electron-electron interaction within a single¬ 
particle framework. Defect states are identified as the 
one-electron Kohn-Sham eigenenergies.® 19 While this 
approach does provide a qualitative picture of the defect- 
derived energy levels, there are two fundamental draw¬ 
backs. The first one is that the one-electron levels of DFT 
have no physical meaning as they relate to the energy 
levels of a fictitious set of non-interacting electrons! 11 ' 
The second issue relates to the band gap problem of tra¬ 
ditional dftP it is well-known that the band gap of 
insulators are significantly underestimated, which place 
uncertainties on the energetic location of the defect en¬ 
ergy levels within the band gap as well as the position of 
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the conduction and valence band edges.® 22 These fac¬ 
tors clearly call into question the emission energies (i.e., 
the energy difference between defect levels and the ap¬ 
propriate band edges) derived from such calculations. 

In the present work, we go beyond the single-particle 
picture of defect levels in PE. Defect levels are charac¬ 
terized in terms of thermodynamic and adiabatic charge 
transition levels that involve total energy differences 
of charged and neutral defect calculations. It is well- 
known that total energy differences provide a formally 
correct and quantitatively better description of energy 
level differences in molecules and solids (including defect 
level placements) IMUiUU Exchange-correlation function¬ 
als more advanced than the (semi)local ones used in the 
past are considered in the present work to critically as¬ 
sess the role of the level of theory adopted. Furthermore, 
many-body perturbation theory computations are under¬ 
taken to accurately determine the band gap of PE and 
its conduction and valence band edge positions. 

The combination of computations undertaken here, in¬ 
cluding the determination and usage of charge transition 
levels (rather than the one-electron levels) and the accu¬ 
rate band edge and gap descriptions allow us to directly 
and quantitatively connect with available luminescence 
data for PE. The present work thus leads to a clear and 
compelling picture of defect states in PE, and provides a 
systematic and reliable pathway for determining the lu¬ 
minescence signatures of various defects in PE as well as 
other insulators. 


II. MODELS AND METHODS 

A. Models 

Perfect crystalline PE consists of chains of singly- 
bonded —CH 2 — groups, of which all the carbon atoms 
display sp 3 hybridization. In this work, 1x2x3 super¬ 
cells of crystalline PE, containing of 24 — CH 2 — groups, 
are constructed with and without defects as the models 
for our calculations. Six kinds of defects, illustrated in 
Fig-i are typically considered to exist in PE chains. A 
carbonyl defect (C=0) is created when the two H atoms 
of a — CH 2 — group are replaced by an O atom doubly 
bonded to the central C atom, which now displays sp 2 
hybridization. If the two H atoms of -CH 2 - are replaced 
by a =CH 2 group, we have a vinyl defect. When each of 
two adjacent C atoms loses one H atom, a double bond 
(-CH=CH-) may be formed between them. The com¬ 
bination of a double bond and a carbonyl is dienone (- 
CH=CH-C=0). In the case of a hydroxyl defect, a H 
atom is replaced by an —OH group. Along the PE chains, 
an alternating pattern of single and double bonds leads 
to conjugated double bonds (-CH=CH-CH=CH-). 


B. DFT computations 

Our DFT calculations were performed with the Vi¬ 
enna ab-initio simulation package (VASP)P4HHD The 
(folded) Brillouin zone of the supercell is sampled by a 
2x2x2 k-point Monkhorst-Pack mestP^ while plane 
waves with kinetic energy up to 400 eV were included 
in the basis set. van der Waals interactions are con¬ 
sidered within the scheme developed by Tkatchenko 
and Scheffler (TS)! 28 -I Relaxed geometries were ob¬ 
tained with Perdew-Burke-Ernzerhof (PBE) exchange- 
correlation (XC) functional,®I which were then used 
to compute the thermodynamic and optical transition 
levels with the PBE and the Heyd-Scuseria-Ernzerhof 
(HSE)|MIi]xC functionals, and the self-consistent many- 
body perturbation theory GW method. In the case of the 
HSE XC functional, we used the HSE06 version with the 
mixing coefficient a = 0.25 and the screening parameter 
uj = 0.2 A -1 . The band gap (E g ) of perfect PE calcu¬ 
lated with PBE, HSE06 and GW methods are 7.00 eV, 
8.39 eV, and 8.95 eV, respectively, with the GW result 
being closest to the experimental value (8.8 eV)P^ 


C. Defect levels theory 

Defect levels are probed by enforcing electronic tran¬ 
sitions, such as those involving a transfer of an electron 
from (to) a defect state to (from) the conduction band 
minimum (the valence band maximum) of PE. Such tran¬ 
sitions can occur in two ways, depending on the time 
scale of the measurement process. Thermodynamic tran¬ 
sitions (involving long time scales) will involve initial and 
final charge states at the respective equilibrium geome¬ 
tries. Adiabatic (or optical) transitions will occur at 
much shorter time scales, and will involve the final charge 
state at the equilibrium geometry of the initial charge 
state. 

A formally correct approach to determine defect lev¬ 
els is via total energy differenc es of PE with defects at 
initial and final charge states.^ 13 33 The local equilib¬ 
rium atomic configurations are different for the defects 
in different charge states,^ as illustrated in Fig. 1(a)- 
A configuration-coordinate diagram of a defect in two 
charge states: 0 (q) and —1 (q') is shown in Fig. [ 2 ] 
(a), in which the minima of the potential energy sur¬ 
faces represent the ground state of the defect in charge 
state q — 0 and q' = — 1. In the case of thermodynamic 
transition, the defect has enough time to relax from its 
initial ground state at q = 0 into its new ground states 
at q' = — 1 or vice versa, as shown using red double ar¬ 
rowhead line in Fig. [ 2 ] (a). In order to get this charge 
transition level, the formation energies Eq(R q ) as a func¬ 
tion of Fermi energy (Ep) for this defect at charge states 
q — 0 and q' = — 1 are computed and shown in Fig. [ 2 ] (b), 
where R q represents the equilibrium structures in charge 
state q. Based on Fig. [ 2 ] (a), the thermodynamic tran¬ 
sition level /i( 0/— 1) corresponds to the crossover point 
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FIG. 2. (Color online) (a) A configuration coordinate dia¬ 
gram for a defect at charge states — 1 and 0 which correspond 
to two equilibrium configurations R q '=-1 and R q = o, shown as 
the minima of the potential energy surfaces (cyan and olive 
lines). The instant transition from 0 state to —1 state gains 
the emission energy E e while the configuration R q =o remains 
unchanged. After a certain amount of time, the —1 state 
evolves to R q >=- 1 , gaining E r 12 from the relaxation process. 
Similarly, the transition from the —1 state to the 0 state by 
absorbing E a happens instantly before the system can relax, 
gaining E r n. (b) Formation energy E q (R q ) as a function of 
Fermi energy ( Ef ) for a defect that has two stable charge 
states: 0 and —1. E q (R q /) refers to the formation energy of 
the defect in charge state q with the equilibrium configura¬ 
tion of charge state q ( R q /), where q and q' are 0 or — 1. 
Solid lines are the formation energies of the relaxed defects 
in each charge state. While the formation energies for the 
defect with frozen atomic configuration of initial charge state 
are presented with dashed lines. Both thermodynamic tran¬ 
sition level fi{ 0/—1) and the relevant optical transition levels 
(M°ot-i) and Mf-Go)) are shown. 


between E i 0 {R q= o) and E^ = _ 1 (R q ' = - 1 ) (solid lines in 
Fig. i (b)). In general, the thermodynamic transition 
levels /i(g/g') is the Fermi energy at which defects in two 
different charge states q and q' are at thermodynamic 
equilibrium and is given by 12 13 33 

q-q' 

Here, E q (R q ) is the formation energy of the g-charged 
defect at its equilibrium structure R q , which can be ob¬ 
tained from DFT calculations. The Fermi energy is taken 
from valence band maximum (VBM) to conduction band 
minimum (CBM) of defect-free PE. In this work, all en¬ 
ergies are referenced to the averaged electrostatic poten¬ 
tial of PE far away from the defect location. Due to 
the unavailability of the total energy with GW method, 
the GW defect levels are corrected by including a term 
AE/vbm = VBMgw — VBM H seo 6 for the VBM shift of 
defect-free PE, expressed as 

{e(q/q'))gw = (m(#/V))hseo6 + AEVbm- ( 2 ) 
The subscripts (HSE06 and GW) are refer to the method 
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FIG. 3. (Color online) Side view of the relaxed PE model with 
defects at different charge states [denoted by (defect) 9 , q — 0, 
— 1 and +1], optimized with PBE functional. The red, white 
and grey spheres are, respectively, O, H and C atoms. The 
bond angles (Z) and bond lengths ( l , given in A) are shown 
in red and black, respectively. 


used to calculate the relevant quantities. 

In the case of optical transition, the atomic configu¬ 
ration of the defect at the initial charge state q is fixed 
even though charge transition exists, i.e., it is a vertical 
transition which is shown using a blue arrowhead line in 
Fig. [ 2 ] (a). Because the optical transition depends on the 
direction of charge transfer, two kinds of optical energies 
are possible, as shown in Fig. [ 2 ] (a): emission energy 
(E e ) for charge transfer from initial state 0 to final state 
— 1 and absorption energy (E a ) for the reverse process. 
The method to determine the optical transition level is 
similar to the previous case, but the energy of the final 
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state is computed using the equilibrium structure of the 
initial state, such as and in Fig. [^(b). 

Therefore, the optical transition level from initial state q 
to final state q' is defined as 


opt . 4 (^)-^(^) 

(y _ q) 


(3) 


where Eq,(R q ) is the defect formation energy in the 
charge state q' corresponding to the equilibrium struc¬ 
ture R q of the initial charge state q. The corresponding 
GW optical transition level is then expressed as 

(M°qXg')^ GW = (ZAZ HSE06 + ^VBM- (4) 


III. RESULTS AND DISCUSSIONS 
A. Defect geometries 

The relaxed structures of PE with defects at differ¬ 
ent charge states q [denoted by (defect) 9 , q = 0,-1 and 
+1] are shown in Fig. [3j The variation of geometry pa¬ 
rameters (bond angle Z and bond length l) of (defect) 9 
indeed shows that the equilibrium structures correspond¬ 
ing to different charge states are different. For example, 
the C=0 bond length of carbonyl increases from 1.228 A 
for the neutral state to 1.315 A for the —1 charge state, 
because the extra electron is localized at the C=0 bond 
and leads to bond weakening. Whereas, compared to 
the (carbonyl) 0 , the decrease of the C=0 bond length in 
(carbonyl) +1 is due to the loss of an electron. As another 
example, we note that a rotation of about 90 degrees can 
be observed for C=C part of (double bond) -1 and (dou¬ 
ble bond) +1 compared to the neutral case. This is be¬ 
cause both positive and negative charges can weaken the 
C=C bond and make the PE chain free to rotate. All 
these relaxations for (defect) 9 directly cause the broad 
luminescence peaks in experiments. Therefore, in this 
work, all levels ranging from which is relevant to 

the emission energy, to n(q/q') are used to compare with 
experimental results. 


B. Charge transition levels 

In order to characterize the defect levels in PE, which 
are responsible for the observed luminescence signatures, 
the thermodynamic and optical transition levels for six 
typical defects in PE have been calculated. The results 
obtained from various methods (solid lines) and avail¬ 
able experimental data are shown in Fig. [4j As can be 
seen, both (0/—1) and (0/+1) charge transition levels are 
obtained for double bond, conjugated double bond and 
vinyl defects. For the carbonyl, dienone and hydroxyl 
defects, just the (0/—1) charge transition level exists, in¬ 
dicating that these three defects at positive charge state 
are not stable within the band gap of PE. The associated 


thermodynamic (i.e., /i(0/—1) and /i(0/+l)) and optical 

/• opt opt opt i opt \ , -J_ • 

(i.e., /i ( o_^_i), M(_ 1 _> 0 ), and M(+i-> 0 )) transitlon 

levels are represented with blue, wine, red, magenta, vi¬ 
olet and orange solid lines, respectively. 

Based on the configuration-coordinate diagram in Fig. 
g(a), we find that the M°oZi) and M°+i are directly 
relevant to luminescence emissions. However, because 
of the broad peak observed in experiments, the levels 
from to fi( 0/—1) or from to /i(0/+1) 

are used to compare with experimental levels, as shown 
by the shaded region in Fig. [4j The experimental levels 
described here, are determined by the difference between 
the CBM and experimental emission energies from Ref. 
H, 34, and [35 To clarify the interpretation, the experi¬ 
mental CBM of defect-free PE is aligned with the GW 
CBM. 

From Fig. [4j we observe that the GW optical transi¬ 
tion levels for carbonyl, dienone, and conjugated 

double bond defects are in good agreement with the ex¬ 
perimental luminescence levels, which suggests that the 
luminescence of such defects originate from a transition 
of an electron from the conduction band to a neutral de¬ 
fect, i.e., (defect) 0 + e (defect) -1 . On the other hand, 
the GW M°+i^o) of vinyl defect is closer to the exper¬ 
imental result compared to This implies that 

the transition of electrons from the conduction band to 
the charge state +1 of vinyl, i.e., (CH = CH 2 ) +1 + e —» 
(CH = CH 2) 0 is most likely the origin of the lumines¬ 
cence center caused by the vinyl defect. However, a small 
discrepancy still remains between computational and ex¬ 
perimental results for the vinyl defect. The reason for 
this disagreement could be because there are uncertain¬ 
ties in the experimental data, or more likely, because the 
present calculations considered defects in a crystalline 
PE environment. The latter aspect may affect the vinyl 
defect the most as this is the bulkiest of defects consid¬ 
ered here. Nevertheless, overall, we are able to achieve 
an acceptable correspondence with available experimen¬ 
tal luminescence data for defects in PE. 

It is also worth noting a few curious details related to 
the computed results. When referenced to the average 
electrostatic potential, the charge transition levels calcu¬ 
lated with PBE agree well with those based on the HSE06 
and the GW method, as shown in Fig. [4] This aspect, 
as already pointed out earlierP 3 reveals that the charge 
transition levels are captured well by the PBE functional, 
if a suitable reference energy is used. 

Unlike the charge transition levels, a remarkable vari¬ 
ation (with the level of theory used) of the band edges 
of PE is observed in Fig. [4j which directly impacts the 
comparison with experimental emission energies (which 
are the energy differences between the defect transition 
levels and the appropriate band edge energies). To show 
the role of band edges, the emission energies of defects 
computed with the PBE, HSE06 and the GW method are 
listed in Table [I| together with available experimental re¬ 
sults. The consistency of GW values and experimental 
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FIG. 4. (Color online) The thermodynamic (i.e., /jl(0/—1) and fi( 0/+1)) and optical (i.e., q\_ 1 ^ 1 ^ 0 )’ ^(oX+i)’ an d 

/i^o)) transition levels are represented with blue, wine, red, magenta, violet and orange solid lines, respectively. The VBM 
and CBM are given with respect to the average electrostatic potential. The possible transition levels which lead to optical 
emissions are represented by the shaded region. The experimental levels are determined by the difference between the CBM 
and experimental emission energies from Ref. El El and [35] For convenience, the experimental CBM of defect-free PE is 
aligned with the GW CBM. For completeness, the Kohn-Sham unoccupied and occupied states are also shown using green and 
cyan dashed lines, respectively. 


TABLE I. The computed emission energies (ranging from CBM — ii°^\ q/ s to CBM — /i(q/q')) and available experimental 
results for defects involved in different charge transfer processes, given in eV. PL, EL, and OL stand for photoluminescence, 
electroluminescence, and other luminescence, respectively. 


Defects 


Charge transfer Optical emission energies 


Experimental values 



Q 



PBE 

HSE06 GW PL EL 

OL 

Carbonyl 

0 


-1 

- 

1.32 - 1.84 2.06 - 3.10 (Ref. 0 2.92 (Ref. Ell 

- 

Dienone 

0 


-1 

1.31 - 1.57 

1.75 - 1.99 2.36 - 3.04 2.25 - 3.54 (Ref. El 2.43 (Ref. EU 

- 

Hydroxyl 

0 


-1 

- 

0.70 

- 

Conjugated dou¬ 
ble bond 

0 


-1 

0.53 - 0.85 

0.91 - 1.23 1.96 - 2.28 2.90 (Ref. EH 

— 

Double bond 

0 


-1 

- 

0.92 - 1.42 

- 

Vinyl 

+ 1 

0 

5.26 - 5.40 

5.91 - 6.05 6.96 - 7.10 - - 4.13 - 

- 6.20 (Ref. 0 


emission energies indicates that an accurate prediction of 
band edges is key to the correct estimation of the emis¬ 
sion energies. 

Based on our results, we propose the following system¬ 
atic procedure to compute emission and absorption en¬ 
ergies accurately: (1) calculation of and n(q/q') 

using the PBE functional; (2) calculation of VBM and 
CBM of pure PE with the GW method; (3) correction 
of and /i(q/q f ) with VBM values from the GW 

method (using Eqs. © and @) by using the average 
electrostatic potential as the energy reference; and fi¬ 


nally (4) the determination of the emission energies as 
the energy differences between appropriate defect transi¬ 
tion levels and band edges. Generally, this procedure can 
be used to predict and assign spectroscopic signatures of 
defects in other organic polymers as well. 


C. Kohn-Sham levels 

To complete the discussion, and to put our work in 
the context of the existing literature pertaining to de- 
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fects and defect states in PE, we also include the one- 
electron Kohn-Sham defect levels for various defects in 
PE in Fig. [4j These are shown using green (unoccupied 
states) and cyan (occupied states) dashed lines in Fig. 
[U While the Kohn-Sham defect levels are in reasonable 
and qualitative agreement with the charge transition lev¬ 
els, some important observations should be made. First, 
even when the PBE and HSE06 average electrostatic po¬ 
tentials are aligned, the one-electron levels do not line up 
(unlike the charge transition levels). This indicates that 
the one-electron levels are strongly dependent on the level 
of theory used. Second, the emission energies that one 
would arrive at purely using the one-electron energies and 
the band edges of the same theoretical treatment would 
be in stark disagreement with experiments. Thus, inter¬ 
preting defect spectroscopic signatures using Kohn-Sham 
levels should be performed with caution. 

IV. SUMMARY 

In summary, we have presented a detailed first- 
principles computational study of chemical defects in 
polyethylene (PE) and have quantitatively determined 
defect-derived emission and absorption energies. A recipe 
for accurately computing and assigning luminescence 
peaks for materials containing defects is proposed. This 
procedure involves going beyond traditional density func¬ 
tional theory calculations within the one-electron frame¬ 
work. It is concluded that defect levels should be deter¬ 
mined in terms of thermodynamic and adiabatic charge 
transition levels (rather than the one-electron Kohn- 
Sham levels) that involve total energy calculations of 
neutral and charged defects, and that band edges should 
be determined accurately using many-body perturbation 
theory. As a consequence of these steps, this work of¬ 
fers a quantitative and satisfactory agreement between 
computed and measured emission/absorption energies for 
polyethylene containing chemical defects. It is suggested 
that this computational scheme may be used to interpret 
the luminescence data of other organic polymers as well. 

V. ACKNOWLEDGEMENT 

This work is supported by a Multidisciplinary Univer¬ 
sity Research Initiative (MURI) grant from the Office of 
Naval Research. 

1 A. Peacock, Handbook of polyethylene: structures: properties, 

and applications (CRC Press, New York, 2000). 

2 R. Patsch, J. Phys. D: Appl. Phys. 23, 1497 (1990). 


3 M. Ieda and Y. Suzuoki, IEEE Electr. Insul. Mag. 13, 10 (1997). 
4 L. Dissado, G. Mazzanti, and G. Montanari, IEEE Trans. Di- 
electr. Electr. Insul. 4, 496 (1997). 

5 D. Ceresoli, E. Tosatti, S. Scandolo, G. Santoro, and S. Serra, 
J. Chem. Phys. 121, 6478 (2004). 

6 G. Mazzanti and G. Montanari, IEEE Trans. Dielectr. Electr. 
Insul. 12 , 876 (2005). 

7 M. Taleb, G. Teyssedre, S. Roy, and C. Laurent, IEEE Trans. 
Dielectr. Electr. Insul. 20, 311 (2013). 

8 N. Allen, J. Homer, J. McKellar, and D. Wood, J. Appl. Polym. 
Sci. 21 , 3147 (1977). 

9 G. Teyssedre, G. Tardieu, and C. Laurent, J. Matter. Sci. 37, 
1599 (2002). 

10 P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964). 

11 W. Kohn and L. Sham, Phys. Rev. 140, A1133 (1965). 

12 C. Freysoldt, B. Grabowski, T. Hickel, J. Neugebauer, G. Kresse, 
A. Janotti, and C. G. Van de Walle, Rev. Mod. Phys. 86, 253 
(2014) 

ia R. Ramprasad, H. Zhu, P. Rinke, and M. Scheffler, Phys. Rev. 
Lett. 108, 066404 (2012) 

14 Q. Yan, A. Janotti, M. Scheffler, and C. G. Van de Walle, Appl. 
Phys. Lett. 105, 111104 (2014). 

15 C. R. Bealing and R. Ramprasad, J. Chem. Phys. 139, 174904 
(2013). 

16 A. Huzayyin, S. Boggs, and R. Ramprasad, IEEE Trans. Di¬ 
electr. Electr. Insul. 17, 926 (2010). 

17 A. Huzayyin, S. Boggs, and R. Ramprasad, IEEE Trans. Di¬ 
electr. Electr. Insul. 17, 920 (2010). 

18 M. Meunier, N. Quirke, and A. Aslanides, J. Chem. Phys. 115, 
2876 (2001). 

19 J. Anta, G. Marcelli, M. Meunier, and N. Quirke, J. Appl. Phys. 
92, 1002 (2002). 

20 A. Seidl, A. Gorling, P. Vogl, J. Majewski, and M. Levy, Phys. 
Rev. B 53, 3764 (1996). 

21 M. S. Hybertsen and S. G. Louie, Rev. Lett. 55, 1418 (1985). 
22 M. Griming, A. Marini, and A. Rubio, J. Chem. Phys. 124, 
154108 (2006). 

23 Q. Yan, A. Janotti, M. Scheffler, and C. G. Van de Walle, Appl. 
Phys. Lett. 100, 142110 (2012). 

24 G. Kresse, Ab initio Molekular Dynamik fiir fliissige Metalle , 
Ph.D. thesis, Technische Universitat Wien (1993). 

25 G. Kresse and Furthmiiller, J. Comput. Mater. Sci. 6, 15 (1996). 
26 G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 (1996). 
27 H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 (1976). 
28 A. Tkatchenko and M. Scheffler, Phys. Rev. Lett. 102 , 073005 
(2009). 

29 J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 
3865 (1996). 

30 J. Heyd, G. E. Scuseria, and M. Ernzerhof, J. Chem. Phys. 124, 
219906 (2006). 

31 A. V. Krukau, O. A. Vydrov, A. F. Izmaylov, and G. E. Scuseria, 
J. Chem. Phys. 125, 224106 (2006). 

32 K. Less and E. Wilson, J. Phys. C 6, 3110 (1973). 

33 C. G. Van de Walle and J. Neugebauer, J. Appl. Phys. 95, 3851 
(2004). 

34 C. Laurent, G. Teyssedre, S. Le Roy, and F. Baudoin, IEEE 
Trans. Dielectr. Electr. Insul. 20, 357 (2013). 

35 T. Ito, D. Kaneko, and Y. Ohki, in 2002 Annual Report Confer¬ 
ence on Electrical Insulation and Dielectric Phenomena (IEEE, 
Quintana Roo, Mexico, 2002) pp. 856-860. 


